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Abstract 

We report on the response of a prototype CMS hadron calorimeter module to 
charged particle beams of pions, muons, and electrons with momenta up to 375 
GeV/c. The data were taken at the H2 and H4 beamlines at CERN in 1995 and 
1996. The prototype sampling calorimeter used copper absorber plates and scintil- 
lator tiles with wavelength shifting fibers for readout. The effects of a magnetic field 
of up to 3 Tesla on the response of the calorimeter to muons, electrons, and pions 
are presented, and the effects of an upstream lead tungstate crystal electromag- 
netic calorimeter on the linearity and energy resolution of the combined calorimet- 
ric system to hadrons are evaluated. The results are compared with Monte Carlo 
simulations and are used to optimize the choice of total absorber depth, sampling 
frequency, and longitudinal readout segmentation. 

PACS Numbers: 29.40.Vj; 29.90.+r; 29.40.Mc; 29.40.Vj Calorimeters 29.90.+r 
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1 The CMS Scintillator Tile-Fiber Hadron Calorimeter 



The hadron calorimeter (HCAL) for the Compact Muon Solenoid (CMS) De- 
tector will be used in the determination of quark, gluon, and neutrino final 
state momenta by measuring the energies and directions of particle jets and 
the missing transverse energy flow. The determination of the missing energy 
flow is crucial in searches for new particles and phenomena, such as possible 
super symmetric partners of quarks and gluons. Adequate granularity, reso- 
lution, and containment of particle showers are essential in attaining these 
performance goals and provide one of the benchmarks in the design of the 
CMS hadron calorimeter. In this communication we report on test beam re- 
sults [1] used in the optimization of the design of HCAL, including choice of 
total absorber depth, sampling frequency, and longitudinal readout segmen- 
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tat ion. 



The central CMS hadron calorimeter [2,3] is located inside the 4 Tesla coil 
of the CMS solenoid magnet (inner diameter of 5.9 m). The central pseudo- 
rapidity range ( \rj\ < 3.0 ) is covered by the barrel and endcap calorimeter 
systems. Segments of the crystal PbW0 4 electromagnetic calorimeter with a 
silicon pre-shower detector (in the endcap region only) are supported by the 
barrel and endcap hadron calorimeters. The combined response of the PbW04 
crystal electromagnetic calorimeter and the hadron calorimeter is used in the 
reconstruction of particle jets and missing energy in the central pseudorapidity 
range. 

The barrel and endcap sampling hadron calorimeters are located inside the 
4 Tesla field of the CMS solenoid. Therefore, the calorimeters are constructed 
with non-magnetic material. The absorber and structural elements are made 
out of cartridge brass (70% Cu/ 30% Zn), and stainless steel plates, respec- 
tively. Cartridge brass is easy to machine and its hadronic interaction length 
is approximately 10% shorter than iron. The active sampling elements are 
3.7 mm thick plastic scintillator tiles with wavelength-shifting (WLS) fiber 
readout. The barrel calorimeter inside the solenoid is relatively thin (about 
7 interaction lengths at rj = 0). To ensure adequate sampling depth for \rj\ 
< 1.4, a hadron outer calorimeter is installed outside the solenoid. The outer 
calorimeter utilizes the solenoid magnet as an additional absorber equal to 
1.4/sin(6 l ) interaction lengths and is used to identify and quantify the contri- 
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bution from late starting showers. 

1.1 Design Requirements 

The design of the central hadron calorimeter requires good hermeticity, good 
transverse granularity, moderate energy resolution, and sufficient depth for 
hadron shower containment. We have chosen a lateral granularity of Ai]xA(f) = 
0.087 x 0.087 for \rj\ < 2.0. This granularity is sufficient to insure good di-jet 
separation and mass resolution. The calorimeter readout is required to have a 
dynamic range from 20 MeV to 3 TeV. The sensitivity at the low end allows 
for the observation of single muons in a calorimeter tower for calibration and 
trigger redundancy. The scale at the high end is set by the maximum energy 
expected to be deposited by a jet in a single calorimeter tower. 

Initial simulations of the CMS calorimeter indicate that a resolution of (Je/E = 
120%/v^E' © 5% for single incident hadrons is sufficient. In this case, the 
energy resolution for a jet of particles between 50 GeV and 3 TeV is not 
degraded by the measurement in the calorimeter, when other fluctuations 
which are inherent in jets [4] are also considered. 

1.2 Magnetic Field Effects 

Particular attention must be paid to the effects of the high magnetic field 
on the response of the calorimeter. The intrinsic shower energy development 
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and the containment of hadron showers have been shown not to be affected 
by the high magnetic field. However, the scintillator itself exhibits an in- 
creased response in a magnetic field. The increased response originates from 
two sources: (1) an intrinsic brightening of the scintillator 5%) for B fields 
above 0.3 Tesla, and (2) a geometric effect due to the increased path lengths 
of soft electrons in the scintillator. The latter effect depends on the specific 
orientation of the magnetic field relative to the calorimeter absorber plates. 

1.3 The Test Beam Program 

A test beam program was initiated in 1994 [5] and was continued in 1995 
and 1996. A test module of moveable brass absorber plates and scintillator 
tile-fiber sampling was exposed to negative hadrons, electrons, and muons 
in the CERN test beam over a large energy range. During this period the 
effects of the following sources on the response of the calorimeter module were 
investigated: 

• Magnetic field (scintillator "brightening" and geometric effects) 

• Absorber thickness (optimization of resolution versus containment) 

• Absorber depth (energy containment) 

• Absorber non-uniformity (presence of the magnet coil) 

• Crystal electromagnetic calorimeter (e/7r effects) 
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Data were taken both with and without a prototype lead tungstate crystal 
electromagnetic calorimeter (ECAL) placed upstream of the HCAL module. 
These data have been used in the optimization of the HCAL design, includ- 
ing the choice of total absorber depth, sampling frequency, and longitudinal 
readout segmentation. 

The test beam data show that the presence of the crystal electromagnetic 
calorimeter in front of HCAL degrades the linearity and resolution of the 
combined calorimeter system. However, a large fraction of this degradation can 
be corrected for by applying constant energy-independent weighting factors 
to the various longitudinal readout segments of HCAL. The response of the 
combined calorimeter system using these optimized weights meets the design 
requirements for resolution and containment of hadron showers. 

The experimental setup, various types of studies, and data sets are described 
in section 2. The measurement of the effects of a magnetic field on the re- 
sponse of HCAL to muons, electrons, and pions is presented in section 3. The 
performance of HCAL and the combined ECAL+HCAL system is discussed 
in sections 4 and 5. Extraction of the intrinsic ratio of the response of HCAL 
to the electromagnetic and hadronic components of the shower, e/h, is pre- 
sented in section 6. The optimization studies leading to the choices for the 
total absorber depth, absorber sampling frequency, and longitudinal readout 
segmentation is presented in section 7. A comparison of the test beam re- 
sults with Monte Carlo simulations is discussed in section 8. Section 9 gives a 
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summary of the results and conclusions. 
2 Experimental Setup 

The combined ECAL+HCAL calorimetric system for CMS has been tested 
in 1995 and 1996 at the H2 and H4 beamlines at CERN. Data were taken 
using beams of muons, electrons, and hadrons, ranging in momenta from 15 
to 375 GeV/c. These test beam prototypes of HCAL are based on the hanging 
file [5] structure, in which copper alloy [6] absorber plates (A/att(Cu)=15.06 
cm) varying in thickness from 2 cm to 10 cm are interspersed with scintillator 
tiles read out with WLS fibers [7]. Each scintillator tile is read out inde- 
pendently with a photomultiplier tube (PMT) located at the end of 10 m 
long optical cable. The total interaction depth of the HCAL prototypes cor- 
responds to 10.3 \ INT (H2, 1995 setup), 8.5 \ INT (H4, 1995 setup), and 10.1 
Xint (H2, 1996 setup). The transverse size of the prototype scintillator tiles is 
64 cmx64 cm. Details of the H2(1995), H4(1995), and H2(1996) experimental 
setups are shown in Figure 1. The longitudinal segmentation of each of the 
modules is described in Tables 1 through 3. 

In the H2(1995) run, the HCAL prototypes have been tested with the detector 
placed inside a large 3 Tesla magnet [8]. The orientation of the magnetic 
field, with B field lines perpendicular to the scintillator planes (magnetic field 
parallel to the beam) corresponds to the endcap configuration of a typical 
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collider detector. The H2(1995) tests include an ECAL module consisting of a 
lead-scintillator sampling calorimeter, with 10 layers of 1.6 cm Pb interspersed 
with 6 mm thick scintillator plates. The transverse size of the lead-scintillator 
sampling ECAL calorimeter for the H2(1995) is 32 cmx32 cm. 

During the H4(1995) run, the combined lead tungstate crystal ECAL+HCAL 
CMS prototype calorimeter has been tested with no magnetic field. The pro- 
totype H4(1995) ECAL detector consists of a matrix of 7x7 PbW0 4 crystals, 
each 23 cm long (25.8 X , 1.1 Xint), and 2 cmx2 cm in transverse dimensions. 
The total transverse size of the ECAL module matrix is 14 cmxl4 cm. The 
crystal ECAL module is placed approximately 50 cm upstream of the front 
face of the HCAL prototype. The PbW0 4 crystals are read out by avalanche 
photodiodes [9] (APD) with a gain approximately equal to 50. 

The combined crystal ECAL+HCAL calorimetric system has been tested 
again in 1996 at the CERN H2 beamline with a magnetic field. The 3 Tesla 
magnet is oriented in such a way that the lines of the B field are parallel to 
the scintillator planes, corresponding to the HCAL barrel configuration. For 
the first 5 interaction lengths of the calorimeter, the magnetic field is uniform 
to within 10%. Figure 2 shows the relative value of the magnetic field, as a 
function of depth in HCAL. 
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2. 1 The Relative Calibration of HCAL Layers 



For the H2(1995) and H4(1995) tests, each scintillator tile is read out by 
20 parallel WLS fibers. Figure 3 shows the details of the design of a single 
scintillator tile. Twenty parallel grooves, spaced every 3.2 cm, are machined 
with a ball-groove tool bit in 4 mm thick SCSN-81 scintillator plates. The 
tiles are painted at the edges using Bicron BC-620 white reflective paint, 
wrapped with white reflective Tyvek [10], and sealed for light leaks with a 
black Tedlar [11]. The WLS fibers are mirrored at one end, placed in the 
grooves, and epoxied in groups of ten into two optical connectors placed at 
the edge of each scintillator tile, as shown in Figure 3. The stainless steel tubes 
are installed inside a 2 mm plastic cover and used to guide a radioactive Cs 137 
source for calibration purposes. The entire package is held together with a set 
of brass rivets. 

For the H2(1996) tests, a different set of scintillator tiles is used. These tiles are 
closer to the tower geometry of the final CMS design. Each tile is segmented 
into a 3x3 matrix read out by nine WLS fibers inserted into a pattern grooves 
and is connected to a single PMT. The scintillator tile layout is shown in Fig- 
ure 4. The tiles are made from 4 mm thick SCSN-81 scintillator produced 
by Kuraray. Multiclad Y-ll WLS fibers, made by Kuraray, 0.83 mm in di- 
ameter, are used. The cross-sectional view of a scintillator tile is shown in 
Figure 5. Scintillator tiles are packaged between 2 mm thick and 1 mm thick 
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plastic covers. The transverse uniformity of response of the scintillator plates, 
measured [12] using a radioactive photon source, has an rms of 4%. 

A schematic view of the CERN H2 test beam setup is shown in Figure 6. 
The calorimeter is placed inside a superconducting magnet, which provides a 
magnetic field of up to 3 Tesla. Each scintillator tile is read out independently 
by a 10 stage, EMI 9902 KA photomultiplier tube. Optical cables carry light 
from the scintillator tiles to the PMT box. In order to avoid large PMT gain 
variations from the fringe fields of the magnet, the PMTs are located approx- 
imately 5 m away from the magnet, behind aim thick iron wall. At this 
location, the fringe field for the magnet at 3 Tesla is approximately 50 Gauss. 
Therefore, each PMT can be adequately shielded by a regular inner //-metal 
shield and an outer, 3 mm thick, soft iron pipe. 

The phototube calibration system is a crucial part of the setup. In order to 
understand the effects of the B field on the phototube gains, several redundant 
calibration schemes to monitor the PMT gain variation are used. The relative 
calibration of the scintillator tiles is established with an accuracy of ~ 3% 
by equalizing the average response of each layer to muons. Figure 7 shows 
an ADC spectrum for 225 GeV/c muons traversing a single scintillator layer. 
The absolute light yield of each scintillator tile, determined using muons, is 
approximately 1.5 photoelectrons per minimum ionizing particle. This light 
yield is sufficiently high such as to not contribute significantly [13] to the 
hadronic energy resolution of HCAL. 
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An independent method of monitoring the gain of each tile is provided by 
a system of radioactive source guide tubes [14]. A radioactive Cs 137 source 
located at the tip of a long stainless steel wire is periodically inserted into 
guide tubes embedded in the plastic cover sheet for each tile. The source is 
moved using a computer controlled source driver. 

The muon calibration and the wire source calibration schemes monitor varia- 
tions in the product of the light yield of the scintillator tiles and the gain of the 
photomultipliers. In order to separately monitor variations in phototube gains 
(from temperature dependence and B fringe field effects), each PMT is also 
connected to a special scintillator block via a set of separate Yll calibration 
fibers. The scintillator block is excited using either a laser beam, or a Cs 137 
source. In addition, the Yll calibration fibers can be excited by a light from a 
blue light emitting diode (LED). The scintillator block is located in a region 
with a low (approximately 50 Gauss) magnetic field, and thus is in a region 
where scintillator brightening from the fringe magnetic field is below 1% [15]. 

3 Calorimeter Performance in a Magnetic Field 

The CMS hadron calorimeter, which uses plastic scintillators, will operate in 
a high magnetic field (4 Tesla). As has been previously observed [15-17], the 
light yield of plastic scintillators placed in high magnetic fields increases by 5 
to 8%. This increase is due to polymer excitation that increases the energy of 
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the short wavelength primary light in the ultra-violet (UV) region. 

In addition, the presence of a magnetic field can affect [5,18] the energy de- 
posited in active parts of a sampling calorimeter. To understand this, the re- 
sponse of HCAL to muons, electrons, and pions has been studied as a function 
of magnetic field strength and orientation for several different sampling config- 
urations. For all magnetic field studies the ECAL module was moved out of the 
beamline. During the H2(1995) tests, the B-field was oriented perpendicular 
to the scintillator planes, corresponding to the endcap HCAL configuration, 
while during the H2(1996) tests, the B-field was parallel to scintillator planes, 
corresponding to the barrel HCAL configuration. 

3.1 B field perpendicular to the scintillator planes (endcap configuration) 

When the magnetic field lines are perpendicular to the scintillator planes, the 
intrinsic light yield of the scintillator is higher than in the case where there 
is no magnetic field. This scintillator brightening effect leads to the same 
overall increased response of the calorimeter to muons, electrons, and pions. 
Figure 8 shows the response of HCAL to pions and electrons as a function of 
B field relative to B=0 Tesla. The observed ~ 5% increase in the response of 
the calorimeter to particles is consistent with the increased light yield of the 
scintillator plastic, as measured by the calibration system using the radioactive 
Cs 137 photon source. 
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Figure 9 shows a comparison of the average shower profile for 300 GeV/c pions 
for two magnetic field settings: B=0 Tesla and B= 3 Tesla. The average energy 
deposition per sampling layer, in minimum ionizing particles (mips), is plotted 
as a function of HCAL absorber depth. Since the gain of each scintillator layer 
is established by equalizing its response to muons independently for the T 
and 3 T data sets, the effect of scintillator brightening is already removed. 
As can be seen in the figure, the shape of the pion shower profiles remains 
unchanged in the presence of a 3 Tesla magnetic field perpendicular to the 
scintillator planes. 

In summary, we conclude that a B field perpendicular to the scintillator planes 
results in an overall increase in the scintillator light yield. This scintillator 
brightening effect can be effectively measured by the calibration system using 
a radioactive 7 source. 

3.2 B field parallel to the scintillator planes (barrel configuration) 

In contrast, when the magnetic field lines are parallel to the scintillator planes, 
the observed average hadron shower profiles are altered. A comparison of the 
average observed shower profiles for 300/c GeV pions as a function of the 
HCAL absorber depth, for B = and B = 3 Tesla, is shown in Figure 10. 
Energy deposited by hadron showers in mips shows an additional increased 
response, in the B = 3 Tesla field, relative to the B= Tesla case. 
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The relative average response of HCAL to 100 GeV/c electrons and pions as 
a function of magnetic field is shown in Figure 11. As indicated by the data, 
the HCAL response to electrons in mips increases by approximately 20% at 
B=3 Tesla, relative to B=0 Tesla. However, the response for pions in mips 
increases by only m 8%. The results are compared to simulations using the 
GEANT [19] Monte Carlo program (10 GeV/c electrons and 50 GeV/c pions). 
As shown in the figure, the simulations reproduce data well. 

Studies using a GEANT Monte Carlo simulation indicate that the increased 
response is due to a geometric effect resulting from a change in the path length 
of low energy electrons (between 1 and 10 MeV) exiting the absorber plates 
and traversing a circular path in the scintillator layer (because of the strong 
magnetic field). The radius of curvature of a 1 MeV/c electron in a 3 Tesla 
magnetic field is approximately 1 mm. The size of this geometric effect is 
expected to be proportional to the strength of the B field and depends on the 
detailed structure [18] and composition of the absorber/scintillator package. 

A GEANT study models the calorimeter response to 10 GeV electrons as a 
function of the air gap between the scintillator package and the upstream 
absorber (see Figure 12) in the presence of a strong magnetic field. The re- 
sults are summarized in Figure 13. The increase in the electron response is 
reduced as the distance between scintillator package and upstream absorber 
is increased. 
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The results presented so far have been taken with the scintillator package 
placed in the absorber gap in the following way: 1 mm plastic + 4 mm scintil- 
lator + 2 mm plastic, with respect to the beam direction. In this configuration 
(configuration A), the scintillator is separated from the upstream- most ab- 
sorber plate by 1 mm ± 0.5 mm of air and 1 mm of plastic. The Monte Carlo 
studies indicated that the effect of a parallel magnetic field depends on the 
distance between the scintillator and the most upstream absorber. Therefore, 
we have investigated the calorimeter response to pions, electrons, and muons 
with the reversed orientation of the scintillator package in the absorber gap 
(configuration B). In configuration B, the position of the scintillator package 
relative to the beam direction is: 2 mm plastic + 4 mm scintillator + 1 mm 
plastic. In the case of configuration B, the scintillator is separated from the 
upstream most absorber plate by 1 mm ± 0.5 mm of air and 2 mm of plastic. 

Figure 14 shows the average electron and pion response of HCAL in mips 
as a function of B field, relative to the response at B= Tesla, measured 
with the configuration B scintillator package. Here, we measure the response 
of electrons at B = 3 Tesla to be only 14 ± 1% higher (relative to B =0 
Tesla), compared to 20% ±1% effect for configuration A. GEANT simulation 
reproduces the electron data well. The additional 1 mm plastic upstream of 
the scintillator in configuration B (vs configuration A) helps range out more 
low energy electrons coming out of the upstream absorber plate. 

In order to minimize the effect of the magnetic field on the response of the CMS 
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hadron barrel calorimeter, one would like to maximize the distance between 
the scintillator and the upstream absorber. Therefore the configuration B of 
scintillator placement in the gaps between absorber plates has been chosen for 
the CMS HCAL design. A system of thin brass (Venetian Blind type) springs 
pushes the scintillator tiles radially outwards, such that the scintillator is 
always positioned closest to the downstream absorber plate. The use of springs 
ensures that the distance between inner absorber and the scintillator plate will 
consist of 2 mm ±0.5 mm of air and 2 mm of plastic. 

Figure 15 shows the ratio of HCAL response to pions as a function of pion en- 
ergy, at 3 Tesla and 1.5 Tesla, relative to its response at Tesla for the two con- 
figurations of the scintillator package inside the absorber gap. The lines shown 
on the Figure 15 indicate the values for 7r(3T)/7r(0T) and 7r(1.5T)/7r(0T) as 
a function of particle momentum. These functions were calculated using the 
Wigmans [21] and Groom [25] parameterizationsQ of the average fraction of 
the electromagnetic component in pion showers. The average fraction of the 
electromagnetic component in pion showers increases as a function of pion 
energy. With e/h of HCAL greater than unity, this increase results in non- 
linearity of HCAL response to pions. For example, the relative response of the 
calorimeter to 300 GeV/c pions is higher by ~ 5%, with respect to 50 GeV/c 
pions. However, a 3 Tesla magnetic field increases the e/h ratio for HCAL by 

1 Detailed study of the average fraction of the electromagnetic component in pion 
showers, F(ttq), and extraction of e/h is presented in section 6. 
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20 % for configuration A and 14% for configuration B, see Figs. 11 and 14. 
Therefore, a 3 Tesla magnetic field is expected to increase the non-linearity of 
HCAL's response to pions between 50 GeV/c and 300 GeV/c by an additional 
« 2-3%. 

The above estimate is relevant only for the subset of single pions which interact 
in HCAL (minimum ionizing in ECAL). However, in the actual CMS detector, 
the common situation is that of a jet of particles depositing a large fraction of 
energy in the ECAL calorimeter. Since the ECAL calorimeter is not sensitive 
to magnetic field effects, the additional non-linearity from the 3 Tesla B field 
is greatly reduced. 

In summary, a B field parallel to the scintillator planes results in an increased 
response due to two effects: (a) the scintillator brightening effect and (b) 
an additional increase in the response to the electromagnetic component of 
hadron showers (geometrical effect). The latter effect is strongly dependent 
on the placement of the scintillator packages with respect to the absorber 
plates of the calorimeter. 

The results reported in this paper are consistent with earlier reports on the 
behavior of plastic scintillators in magnetic fields [15,17] and on the influ- 
ence of magnetic fields on response of scintillator based calorimeter [16]. Note 
however, that in the latter publication the conclusion that electron response 
follows the dependence of the light yield ("brightening effect") is true only for 
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magnetic fields below 0.3 Tesla. As shown in this report, in the presense of par- 
allel magnetic fields above 0.3 Tesla, the geometric effect leads to additional 
increase in response of the sampling calorimeters to electrons. 

We also note that the origin of the latter effect is due to the relative geometry 
of the absorber with respect to the active elements of the calorimeter and 
is not sensitive to any special property of the scintillator as the active sam- 
pling medium. Therefore, similar effects are expected for any type of sampling 
calorimeters, with similar magnetic field/absorber configurations, independent 
of the choice of material for the active sampling medium. 

4 Performance the HCAL calorimeter with a Pb-Scintillator ECAL 

During the H2(1995) running period, a Pb-scintillator ECAL calorimeter was 
placed in front of the HCAL calorimeter. Linearity of the response and energy 
resolution of the combined system was studied for electrons, pions, and muons 
with momenta from 50 to 300 GeV/c. 

In the following discussion we consider two different samples of pion data. The 
first, which we refer to as " mip-in-ECAL" pions, includes only pions which are 
minimum ionizing in ECAL, i.e. do not interact in ECAL and begin interacting 
in HCAL. This sample is selected by requiring that the energy deposition in 
ECAL is consistent with minimum ionizing deposition, i.e. is less than 2.5 
GeV. The second sample, which we refer to as the "full pion sample", includes 
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all pions, i.e. interacting in either ECAL or HCAL. 

The absolute energy scale of the HCAL calorimeter is set using 300 GeV/c 
" mip-in-ECAL" pions. The average pulse height deposited by muons in HCAL 
corresponds to approximately 3 GeV of equivalent hadronic energy. Figure 16 
shows the energy deposited by 225 GeV/c muons in HCAL. A pedestal trigger 
is a random trigger during the beam. This pedestal distribution is shown as 
a dashed line in Fig. 16. The width of the energy distribution for random 
pedestal triggers is equivalent to 80 MeV. Hence, the contribution of electronic 
noise to the energy resolution of HCAL is negligible. 

The absolute energy scale of the ECAL calorimeter is set using 50 GeV/c elec- 
trons. The energy response of the Pb-scintillator sampling ECAL calorimeter 
2.9 X sampling) for electrons with momenta in range 50 to 150 GeV/c is 
linear to within 1%. The fractional energy resolution for 50 and 100 GeV/c 
electrons is 7.2% and 6.2%, respectively. 

The energy response of the combined Pb-scintillator ECAL+HCAL calorime- 
ter system is calculated as follows: 



Here, Escale and H scale are constants which define the absolute energy scales 
of ECAL and HCAL, and ADCi are the ADC values of each readout layer 
of ECAL and HCAL. The relative calibration constants of each layer, /ij, 




/ Escale + ^ 




ADCj * wj 



/Hscale (1) 
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correspond to the average muon response of a layer expressed in units of ADC 
counts. The coefficients Wi are proportional to the arithmetic mean of the 
interaction lengths of the absorber plates upstream and downstream of each 
scintillator layer (the "Simpson" approximation formula), w = 0.5x(t upstream + 

^downstream ) • 

Scatter plots of the energies in HCAL vs ECAL for 50, 100, 150, and 300 
GeV/c particles in the hadron beam are shown in Figure 17. As indicated 
in the figure, the nominal hadron beam, especially at the 50 GeV/c and 100 
GeV/c beam tunes, contains a large fraction of electrons. Pion induced events 
are selected by requiring that energy deposited in HCAL is at least 0.5 GeV. 
In order to remove particles that interact upstream of ECAL, the energy 
deposition in the first ECAL layer is required to be consistent with that of a 
minimum ionizing particle. 

The total (ECAL+HCAL) energy distributions for "mip-in-ECAL" pions are 
shown in Figure 18. The total energy distributions for the "full pion sample" 
are shown in Figure 19. The reconstructed energy distributions are well de- 
scribed by Gaussian fits for both sets of data. Figure 20 shows the linearity 
versus energy for the calorimeter response for "mip-in-ECAL" pions, and for 
the "full pion sample" using the Pb-scintillator ECAL+HCAL calorimetric 
system. The average response of the combined Pb-scintillator ECAL+HCAL 
system to the "full pion sample" is approximately 1% lower than the response 
the calorimeter for the subset of "mip-in-ECAL" pions. 
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A comparison of the relative energy resolutions of the calorimeter for "mip- 
in-ECAL" pions and for the "full pion sample" is shown in Figure 21. The 
energy resolutions for these two cases are comparable. At low energies, the 
energy resolution of the "full pion sample" is slightly improved because of the 
finer sampling of the Pb-scintillator ECAL. 

The fractional energy resolutions are parameterized by the following function: 



Here E is the particle energy in GeV, and the symbol © implies that the 
stochastic and constant terms in the resolution are combined in quadrature. 

The following values for the stochastic term and the constant term are ex- 
tracted from the data: 



a E /E — (stock. term)j\pE © (const, term) 



(2) 



,pions mip in ECAL 



IE = (93.8 ± 0.9)%/VE © (4.4 ± 0.1)% 



(3) 



for "mip-in-ECAL" pions, and 



all pions 



IE = (82.6 ± 0.6)%/ VE © (4.5 ± 0.1)% 



(4) 



for the " full pion sample" . 
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5 Performance of the HCAL calorimeter with a PbWC>4 crystal 
ECAL 

In this section the linearity of response and energy resolution of the HCAL 
calorimeter with a PbW0 4 crystal matrix ECAL are discussed. The perfor- 
mance of the combined PbW0 4 ECAL+HCAL calorimetric system was first 
investigated during 1995 in the H4 beamline. 

The calibration of the ECAL crystals uses 50 GeV/c electrons directed into the 
center of each crystal. Figure 22 shows the energy response of ECAL for 25, 50, 
100, and 150 GeV/c electrons using the 7x7 ECAL matrix sum. The energy 
resolutions of the crystal ECAL at these energies are found to be 0.6 GeV, 0.7 
GeV, 1.0 GeV, and 1.4 GeV, corresponding to fractional energy resolutions of 
2.4%, 1.4%, 1.0%, and 0.9%, respectively. During these tests, the electronic 
noise of the 7x7 crystal matrix energy sum had a rms width equivalent to 
440 MeV. Much better performance of the crystal ECAL calorimeter has been 
attained [20] in other tests which were dedicated to special ECAL studies. 

The energy response of the combined crystal ECAL+HCAL calorimeter sys- 
tem is calculated as follows: 

ECAL a nr HCAL a nr . . 

E TOT = ]T f^l/Escale+ ]T W] /H scale (5) 

i e i j /A? 

Here, Escale and Hscale are constants which define the absolute energy scales 
of ECAL and HCAL, ADCi are the ADC values of each readout layer (tower) 
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of the HCAL (ECAL) calorimeter. The relative calibration constants of each 
tower of ECAL, e^, correspond to the average 50 GeV/c electron response of a 
tower expressed in units of ADC counts. The coefficients Wj and fij are defined 
as before (Eqn. 1). 

Figure 23 shows a scatter plot of the energy deposition in HCAL versus ECAL 
taken during the H4(1995) test beam run with pions with momenta ranging 
between 15 and 375 GeV/c. Unlike in the case with the Pb-scintillator ECAL, 
the data points are not centered on a straight line corresponding to Eecal + 
Ehcal = Pbeam- The energy response of the HCAL calorimeter to "mip-in- 
ECAL" pions is shown in Figure 24. The response of the combined crystal 
ECAL+HCAL calorimeter to the "full pion sample" is shown in Figure 25. 
The distributions for the combined crystal ECAL+HCAL system deviate from 
a Gaussian shape, especially at low energies. 

As demonstrated earlier, the H2(1995) data with a Pb-scintillator ECAL show 
good linearity for the combined ECAL+HCAL calorimeter system. However, 
as can be seen in Figure 26, the H4(1995) data with a PbW0 4 crystal ECAL 
is quite non-linear. In addition, the average response of the combined crystal 
ECAL+HCAL calorimeter to "full pion sample" is approximately 10% lower 
than the response to pions which interact in HCAL only. 

Figure 27 shows the fractional energy resolutions of the calorimeter for the 
"mip-in-ECAL" pions for the "full pion sample". The energy resolutions of 
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the combined crystal ECAL+HCAL calorimeter system are significantly de- 
graded, especially at high energies. These results are attributed to the large 
difference between the response of the crystal ECAL to electrons and hadrons, 
i.e. the large e/h (the ratio of the response to the electromagnetic and hadronic 
components of the shower) of the PbW0 4 crystals. This large e/h results in 
an overall response of the calorimeter which is very sensitive to fluctuations 
in the initial electromagnetic component of hadronic showers. 

A fit to data for the fractional energy resolution yields the following values for 
the stochastic term and constant term parameters: 

al nons mip in ECAL ^ = ^ ± Q ^%^ © (4.0 ± 0.1)% (6) 

for the case of " mip-in-ECAL" pions , and 

a aii pions i E = (127 ± , 7)%/v /£ ( 6 5 ± .i)% ( 7 ) 

for the case of the "full pion sample". Note that the first fit ("mip-in-ECAL" 
data , eqn. 6) is consistent with the H2(1995) "mip-in-ECAL" data (eqn. 3). 

6 Response of HCAL to electrons and the e/h ratio for HCAL 

The linearity and energy resolution of hadron calorimeter depends [21] on 
the intrinsic e/h. The average fraction of the energy in the electromagnetic 
component of hadron induced showers, F(tt°), increases as a function of in- 
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cident energy. For non-compensating calorimeters ( e/h ^1), this implies a 
non-linear hadron energy response. The event-by-event fluctuations in F(tt°) 
contribute to variations in the reconstructed shower energy and dominate the 
energy resolution for hadrons at high energies. 

In order to study the relative response of HCAL to electron and pions, data 
were taken with ECAL removed from the beamline. In order to reject electrons 
present in the hadron tunes from the sample, a minimum ionizing signal was 
required in the first three sampling layers of HCAL. 

The studies of the energy response of HCAL to pions indicate that the HCAL 
calorimeter is somewhat non-linear. There is ~ 9% increase in the relative 
response of HCAL to pions between 15 and 375 GeV/c, as shown in Figure 28. 
The total interaction length [f]of the H4(1995) HCAL module is approximately 
8 \int- At very high energies, an 8 Xint calorimeter is not expected to fully 
contain the hadronic shower. The correction for longitudinal leakage for an 
8 Xint calorimeter is estimated [22] to be 1.5% at 120 GeV/c and 3.5% at 
375 GeV/c, using CCFR [23] and NuTeV [24] test beam data. Therefore, 
the intrinsic non-linearity of HCAL is somewhat larger if the corrections for 
leakage at high energies are applied. 

The response of the H4(1995) HCAL module to electrons is linear at a level 

2 Taking into account the requirement that energy deposition in first three HCAL 
counters is consistent with minimum ionizing particles. 
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of 2% and is approximately 15% higher than the average response to pions as 
shown in Figure 28. The absolute energy scale of HCAL is set (as before) by 
the 375 GeV/c pion data point. 

The electron and pion data can be used to extract e/h by applying the fol- 
lowing expressions: 

7T = F(tt°) x e + (1 - F(tt )) x h (8) 

1 = 1 + (f ~ 1) x F(vr°) (g) 

e h 

Here e is the response of the calorimeter to the electromagnetic components 
of hadron showers and h is the response to the hadronic component of pion 
showers, n is the response to pions and F(7r°) is the fraction of the energy in 
electromagnetic component of pion showers. 

The extracted value of e/h depends on the assumed parameterization for F(7r°) 
as a function of energy. The Wigmans [21] formula for F(tt°) (eqn. 10) increases 
with energy (ss log(E)) and becomes non-physical at very high energies. The 
Groom parameterization (eqn. 11) is extracted from a CALOR Monte Carlo 
simulation [25] and uses a power law dependence. 

F w (tt°) = 0.11 x ln{E) (10) 

F g (tt°) = 1 - (E/0.96) - 816 - 1 (11) 

Using the above two parameterizations, we extract values of e/h of 1.41 ± 
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0.01 (Wigmans) and 1.51 ± 0.01 (Groom). The best fit to our data is shown 
in Figure 29. 

The above results can be compared with similar analyses for other iron- 
scintillator sampling calorimeters. Extraction of e/h for an iron-scintillator 
sampling calorimeter has also been done [26] by CDF. The CDF endcap hadron 
calorimeter consists of 22 layers of 5 cm iron absorber plates and 6 mm scin- 
tillators. The CDF results yield e/h 1.34 ± 0.01 (Wigmans) and 1.42 ± 0.015 
(Groom). ATLAS [27] measures e/h= 1.34 ± 0.03 assuming the Wigmans 
parameterization for F(ir°). The ATLAS hadron calorimeter uses iron as an 
absorber with the tiles are oriented in the radial direction. 

The NuTeV Collaboration reports a e/h ratio for 10 cm iron sampling calorime- 
ter [24] of 1.08 ± 0.01 using the Groom parameterization. Note that e/h is sen- 
sitive to the sampling fraction and other geometrical effects. The CCFR/NuTeV 
calorimeter uses 2.5 cm thick liquid scintillator counters, which are clad with 
acrylic and water bags. The cladding and the thick scintillators in CCFR/NuTeV 
greatly reduce the response to electrons relative to the response to pions. Also, 
the water and thick scintillators tend to increases the rate of conversion of low 
energy neutrons to protons in the active medium. 
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7 Optimization of the design of HCAL 

The conceptual design of the barrel HCAL was a 5.3 Xint thick calorime- 
ter. The inner half of the calorimeter had 3 cm sampling (first readout seg- 
ment, HI) while the second readout segment consisted of the outer half of the 
calorimeter with 6 cm sampling. 

For the H2(1996) tests, the prototype HCAL module was segmented into 27 
readout layers. Using these data we simulated various sampling configurations 
and studied the performance of HCAL as a function of total interaction length 
and sampling frequency. These data allowed us to optimize the calorimeter re- 
sponse to pions and jets, while taking into account the existing design (e.g. 
geometrical and readout) constraints imposed on HCAL within the CMS de- 
tector. 

Important design choices that were made based on these data are: 

• the absorber sampling thickness; 

• the depth of the HCAL inside the magnet (in interaction lengths); 

• the longitudinal segmentation of the inner HCAL into HI and H2 readout 
segments; 

• possibility of adding scintillator layers outside of the magnet to create an 
outer calorimeter. 
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7.1 Choice of the longitudinal segmentation for HCAL 

Studies performed prior to the H2(1996) beam tests resulted in no compelling 
argument to set the optimal partition between the two readout segments of 
inner HCAL. However, the 1996 test beam data indicate that the proper par- 
tition is that which is most useful in correcting for large e/h response of 
the ECAL crystals. To compensate for a large e/h of the combined ECAL 
+ HCAL system, we chose to make a novel longitudinal segmentation, with 
the first depth segment, HI, reading out only the first scintillator layer. The 
remaining scintillator layers inside the magnet are combined into the second 
readout, H2. The reason for this choice is the following. A large HI signal 
indicates that a significant amount of hadronic energy has been deposited in 
ECAL and is underestimated because of the large e/h of the crystals. The in- 
formation from HI can be used to make an effective correction for the hadronic 
non-linearity of the ECAL crystals. 

7.2 Response of the combined crystal ECAL+HCAL calorimeters 

Figure 31 and 32 show the energy response and fractional energy resolutions 
of the combined PbW04 crystal ECAL+HCAL calorimeters to pions. The 
inner HCAL(5.3 Xint) consists of two independent readouts: HI (following 
a 2 cm Cu plate) and H2 (thirteen 6 cm Cu samplings). A single readout 
outer calorimeter (HO) consists of three samplings: (1) the first sample is 
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immediately after the magnetic coil (which is mimicked by 18 cm of Cu in the 
test beam ),(2) a second 22 cm Cu sampling layer, and (3) a third 16 cm Cu 
sampling layer. The combined coil and HO samplings correspond to a total of 
2.5A/ATT- 

We have investigated two possible approaches to correct for the degradation 
of the performance of the combined crystal ECAL+HCAL calorimeters. Both 
of these methods make use of the segmented readout of HCAL. The first 
approach, called passive weighting, increases the weight (a) of the first (HI) 
HCAL readout segment, where a is an energy independent constant. Note 
that a =1 corresponds to using the standard "absorber Simpson's weighting" 
used elsewhere in the calorimeter. 

Etot = Eecal + ol x E m + E H2 + E HO (12) 

The overall linearity and fractional energy resolution of the combined crystal 
ECAL+HCAL system for 300 GeV/c pions as a function of the parameter 
a are shown in Figure 30. Both the linearity and energy resolution of the 
combined crystal ECAL+HCAL system are improved for the value of a — 
1.4. 

The second approach, called dynamic weighting is an event-by-event correc- 
tion. It depends on the fraction of the energy deposited in the first readout 
segment of HCAL immediately downstream of the crystal ECAL. The dy- 
namic weighting effectively allows one to have an energy dependent correction 
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(for single particles) for the low response to pions which interact in ECAL. 

Etot = (1 + 2 x f(Hl)) x E EC al + E m + E H2 + E HO , (13) 
f(Hl) = E(H1)/(E(H1) + E(H2) + E(HO)), f(Hl) < 0.1 (14) 

With either the passive or dynamic weighting, the nonlinearity and resolu- 
tions for pions with energies between 30 and 300 GeV/c are improved. With 
the passive weighing, the fractional energy resolution of of the combined 
ECAL+HCAL calorimetric system can be described by the function ge/E 
= 122%/ ' \[E © 5%. Note that while the passive weighting can be applied to 
single particles and jets, the dynamic weighting may introduce high energy 
tails in the case of particle jets. 

Techniques to overcome the problem of the different non-compensation proper- 
ties (i.e. e/h) of the electromagnetic and hadronic compartments for combined 
calorimetry systems have been studied [28] by the ATLAS calorimetry group. 
The algorithms proposed by the ATLAS group depend on the energy frac- 
tion deposited in the electromagnetic compartment and therefore can be only 
applied to single particle energy reconstruction. 

7. 3 Studies of the total HCAL absorber depth 

We use the 27 sampling layers of the H2(1996) test beam module to study 
the calorimeter performance as a function of total depth. Figure 33 shows the 
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average pion shower profiles as a function of total absorber depth for 50, 100, 
150, and 300 GeV/c pions. As shown in the figure, the average pion shower 
profiles extend significantly beyond 7 Xint, in particular at high pion energies. 
Fluctuations in the leakage for high energy pions also become large. 

To collect this leakage energy, the barrel HCAL has been augmented with an 
additional outer calorimeter (HO) consisting of a single layer of scintillators 
beyond the solenoid magnet (the solenoid thickness is 1.4 Xint)- For pseudo- 
rapidity rj less than 0.4, a layer of iron (thickness = 18 cm = 1.1 Xint) has 
been instrumented with an additional scintillator layer. For either case (one 
or two layers of scintillator), the scintillators cover the same solid angle as the 
interior calorimeter and are read out as a single depth segment. The combined 
solenoid + iron of the HO corresponds to an additional 2.5 Xint- 

The fraction of 300 GeV/c pions with reconstructed energy less than 200 GeV 
(approximately 3 o below the mean, or 100 GeV of missing energy) is shown in 
Figure 34. The four points correspond to: 5.9 Xint, the HCAL alone; 7.0 Xint, 
HCAL + ECAL; 9.5 Xint, HCAL + ECAL + HO; and 11 X mT which is 
the total thickness of the 27 layer test beam module. We see that for the 
CMS barrel design, at 9.5 Xint, less than 2% of the pions are catastrophically 
mismeasured. To increase the interaction length of the calorimeter to 9.5 Xint 
we added two additional absorber plates inside the magnet with respect to the 
conceptual design. This was achieved by reducing the inner radius of barrel 
HCAL. 
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7.4 Optimization of Absorber Sampling Thickness 



The fractional energy resolution has been investigated for the following three 
choices for the inner HCAL absorber samplings: a) 3 cm Cu sampling for the 
first eight layers followed by 6 cm Cu sampling, b) 6 cm Cu uniform sampling, 
and c) 12 cm Cu uniform sampling. The data, shown in Figure 35 indicate 
that the energy resolution is not dominated by the sampling fluctuations in 
HCAL. A factor of two decrease in the sampling frequency ( 3 cm/6 cm to 
6 cm uniform sampling) for HCAL does not result in a noticeable degradation 
of the energy resolution of the combined detector system. In the case of 12 cm 
Cu uniform sampling, the degradation in energy resolution is noticeable, but 
does not scale with y/i, where t is the thickness of the absorber plates. For 
the final design, we chose uniform 5 cm sampling which does not degrade the 
energy resolution and puts more absorber inside the magnet, relative to the 
conceptual design. 



7.5 Final Design of CMS Barrel HCAL 



The final design of the CMS barrel HCAL consists of 17 absorber/scintillator 
samples with 5 cm absorber thickness. The total thickness of absorber at 
90 degrees is 5.9 Xint- With the ECAL included, the total thickness inside 
the magnet is 7.0 Xint- Exterior to the magnet, there is a 2.5 Xint outer 
calorimeter. Each projective tower has 3 readouts in depth: the very first 
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scintillator layer (HI); the remaining 16 scintillator layers inside the magnet 
(H2); and the outer calorimeter (HO). 

8 Monte Carlo simulation of the test beam results 

Several hadron shower generators are available within the GEANT [19] frame- 
work. These include GHEISHA [29] , GFLUKA (which is an implementation of 
FLUKA [30] within GEANT), and GCALOR (which is an implementation of 
CALOR [31] within GEANT). Although GHEISHA is native to the GEANT 
program, GFLUKA and GCALOR are imperfect implementations of the orig- 
inal FLUKA and CALOR programs. Both of these programs have been known 
to produce somewhat different results [32] than the original generators. 

GEANT is used in various studies for evaluation of the calorimeter design for 
CMS. In order to verify those simulations and to understand their limitations, 
GCALOR is used to simulate the H2(1996) test beam data. It also serves 
as a reference for comparison with the other generators. Our aim is first to 
anchor the Monte Carlo model to the ensemble of test beam data. Once it is 
so constrained, it is assumed that it can be used to make small extrapolations 
to model the final CMS calorimeter system. 

Details of the crystal ECAL and HCAL test beam geometry are implemented 
in the GCALOR simulation. For ECAL, this includes a 7 x 7 matrix of indi- 
vidual crystals surrounded by copper blocks as well as mechanical and cooling 
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structures. The HCAL geometry includes the layer structure of copper plates, 
scintillator, plastic cover plates, and air gaps on both sides of each scintil- 
lator. In order to take into account all experimental effects, the transverse 
beam profiles are simulated using information from the test beam tracking 
chamber. Electronic noise and photo-statistics effects are simulated based on 
the measured distributions of pedestals, electron, and muon signals in ECAL 
and HCAL. The longitudinal light collection efficiencies in the central nine 
crystals are included in the simulation. The energy cut values in the GEANT 
simulation are set to the GEANT default values, 1 MeV for electrons and 10 
MeV for hadrons. 

In order to compare response of pions interacting in HCAL, we have removed 
ECAL from the beam. For this dataset, the test beam data are in good agree- 
ment with the results of the GCALOR Monte Carlo simulations. Good agree- 
ment is observed in the longitudinal shower profile (Fig. 36), the pion response 
versus energy (Fig.37), and in the fractional pion energy resolution (Fig. 38). 

A comparison of the linearity and fractional energy resolutions for the "full 
pion sample" using the combined crystal ECAL+HCAL system for test beam 
data and the GEANT simulations are shown Figures 39 and 40. The GEANT 
simulations, which include all experimental detector effects such as electronics 
noise, are in good agreement with the test beam data. Also shown in Figure 40 
are the results of a MC simulation of the crystal ECAL+HCAL energy resolu- 
tion excluding test beam detector effects, such as the ECAL electronic noise, 
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and the energy leakage from the small prototype ECAL. As can be seen in 
Figures 40, the detector effects present during the H2(f996) test beam data 
taking significantly degrade the fractional pion energy resolution, especially 
at the low energies. These test beam detector effects are not expected to be 
present in the final CMS configuration. 

Comparisons of the fractional energy resolution simulated by GCALOR and 
other GEANT hadron simulators are shown in Figure 41. The simulations 
based on GHE1SHA predict somewhat worse resolutions, while simulations 
based on GFLUKA+MICAP predict much better resolutions than predicted 
with GCALOR and observed in the data. 

9 Conclusions 

Comprehensive tests of the performance of prototype CMS central hadron 
calorimeters have been done in the H2 and H4 beamlines at CERN. Data 
were taken with both a stand-alone HCAL calorimeter and with HCAL in 
combination with an upstream ECAL calorimeter. 

One of the primary objectives of the HCAL test beam studies was to investi- 
gate the calorimeter performance in the presence of perpendicular and parallel 
magnetic fields. A high magnetic field changes the response of the calorimeter 
in two different ways: (1) the overall light yield of the scintillator is increased 
in a high magnetic field and (2) the field affects the observed energy deposition 
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of electromagnetic components of showers when it is parallel to the scintillator 
plates (i.e. perpendicular to the particle direction). For a collider experiment 
with a solenoid magnet, the magnetic field is parallel to the calorimeter plates 
in the central part of the detector (barrel configuration) and is perpendicular 
to calorimeter plates in the large r\ region ( endcap configuration). 

When the magnetic field is perpendicular to the scintillator planes (endcap 
configuration), only an intrinsic increase of the light yield of the scintillator 
of approximately 5-8% is observed relative to the case with no magnetic field. 
This effect leads to the same overall increased response of the calorimeter to 
muons, electrons, pions, and 7 rays from a radioactive source. Therefore a 
calibration source can be used to track and correct for this effect. 

An additional geometric effect leads to an increased response of the calorimeter 
to the electromagnetic component of showers. This effect occurs in the case 
in which the magnetic field lines are parallel to the scintillator planes (barrel 
configuration) and originates from an increase in the geometrical path length 
of low energy electrons in a magnetic field. Since this effect is not present for 
the calibration source, in situ (B field on) calibration is required for hadron 
barrel calorimeter. 

The size of this effect is approximately proportional to the strength of the B 
field and depends on the detailed structure and composition of the absorber 
and scintillator planes. However, the effect is expected to be present in any 
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sampling calorimeter (independent of readout technology) situated in a mag- 
netic field which is parallel to the readout planes (barrel configuration). 

Due to the non-compensating nature of the lead tungstate crystal ECAL the 
linearity and energy resolution of the combined crystal ECAL+HCAL is worse 
than for HCAL alone. Therefore, improvements in the linearity and resolution 
using weighting of the two longitudinal readouts (HI and H2) of HCAL have 
been investigated. Using a passive weighting, the fractional energy resolution 
of the combined crystal ECAL+HCAL calorimetric system is described by 
the function a E / E—122% / 'y/~E © 5%. This performance is achieved by a novel 
longitudinal segmentation of HCAL with HI being a single layer immediately 
behind ECAL. We conclude that for combined ECAL+HCAL calorimeter sys- 
tems with an ECAL with a large e/h, a very thin first depth segment in HCAL 
can be used to largely correct for the resultant non-linearity and degradation 
in energy resolution. Monte Carlo studies of the CMS detector in a collider en- 
viornment indicate that with the above performance, the energy resolution for 
jets is not dominated by the energy measurement in the hadron calorimeter, 
but by other fluctuations which are inherent in jets [4]. 

We find that the average longitudinal hadron shower profiles extend past the 
inner HCAL located inside the magnetic coil. Therefore, the CMS central 
calorimeter design includes instrumenting an outer calorimeter (HO) outside 
the coil to measure the energy in the tail of high energy hadronic showers. 
The addition of HO leads to a total ECAL+HCAL+HO depth of at least 9.5 
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Xint for almost the entire 77 range spanned by the CMS hadron calorimeters. 

Various GEANT based simulations predict the impact of the performance of 
HCAL on a variety of potential physics searches in the CMS detector. These 
simulations extrapolate the performance of HCAL from the test beam config- 
uration to the final configuration chosen for CMS HCAL design. The Monte 
Carlo programs have been shown to successfully simulate various test beam 
setups, including hadron shower calorimetry in a setting with both crystal and 
copper-scintillator detectors in a strong magnetic field. The best description of 
the data is provided by GEANT with the GCALOR module for the generation 
of hadron showers. 
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Layer number 


Absorber thickness 


Scintillator thickness 


ECAL 1-10 


1.6 cm Pb 


6 mm SCSN-38 


HCAL 1-9 


5 cm Cu 


4 mm SCSN-81 


HCAL 10-20 


10 cm Cu 


4 mm SCSN-81 



Table 1 

Longitudinal segmentation of ECAL and HCAL in the H2(1995) setup. Here a Pb- 
scintillator sampling calorimeter is used for the ECAL section. 
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Layer number 


Absorber thickness 


Scintillator thickness 


T T r \ AT -I 

HCAL 1 


2 cm Cu 


4 mm SCSN-81 


HCAL 2-10 


d cm Cu 


4 mm SCSN-81 


HCAL 11-19 


b cm Cu 


4 mm SCSN-81 


t t /"I At c\r\ 

HCAL 20 


8 cm Cu + 29 cm Al 


4 mm SCSN-81 


T T /"I AT £"\ 1 

HCAL 21 


8 cm Cu 


4 mm SCSN-81 


HCAL 22 


8 cm Cu 


4 mm SCSN-81 


HCAL 23 


10 cm Cu 


4 mm SCSN-81 


HCAL 24 


10 cm Cu 


4 mm SCSN-81 



Table 2 

Longitudinal segmentation of HCAL in the H4(1995) setup. Here the ECAL module 
consists of a 7x7 matrix of PbWC-4 crystals. 
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T 1 

Layer number 


Absorber thickness 


Scintillator thickness 


HCAL 1 


2 cm Cu 


4 mm SCSN-81 


HCAL 2-7 


3 cm Cu 


4 mm SCSN-81 


HCAL 8-21 


6 cm Cu 


4 mm SCSN-81 


HCAL 22-27 


8 cm Cu 


4 mm SCSN-81 



Table 3 

Longitudinal segmentation of HCAL in the H2(1996) setup. Here, the ECAL module 
consists of a 7x7 matrix of PbW04 crystals. 
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CMS HCAL PROTOTYPES 



H2(1995) SETUP: Pb/Sc in. ECAL MODULE 

eca£ 



Transver $p 

Size 
32cm x 3 2cm 



B FIELD 



x 5cm Cu + 11 x 10cm Cu 



H4(1995) SETUP: PbW04 CRYSTAL ECAL MODULE 



ECAir 



50cm ga 



COIL 



NO B FIELD 



2cm Cu+ 9x3cm Cu+ 9x6cm Cu+ 8cmCu+ 29cm A1+ 2x8cm Cu+ 2x1 



H2(1996) SETUP: PbW04 CRYSTAL ECAL MODULE 



5 cm g. 



EC Air 



Transverse 
Size 



f B FIELD 



14cm x 14cm2cm Cu + 6 x 3cm Cu + 14 x 6cm Cu + 6 x 8cm Cu 



Fig. 1. Schematic drawing of the calorimeter modules used during the 1995 and 
1996 tests in the H2 and H4 beamlines. A lead-scintillator sampling ECAL detector 
is used in the H2(1995) setup. A 7x7 matrix of PI3WO4 crystal ECAL is used in 
the H4(1995) and H2(1996) setups. 
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Fig. 2. Relative strength of the magnetic field, as a function of HCAL depth (in 
cm of Cu) for the CMS H2(1996) test beam run. 
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Fig. 3. Design of the 64 cm x 64 cm scintillator tile for the 1995 CMS Test Beam 
HCAL module. 
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Fig. 4. Design of the 64 cm x 64 cm scintillator tile for the 1996 CMS Test Beam 
HCAL module. 
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g. 5. Cross sectional view of the scintillator tile for the CMS Test Beam HCAL 
dule. The dimensions in the drawing are shown in mm. 
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CMS CERN H2 TEST BEAM - TOP VIEW 




Fig. 6. Top view of the CMS test beam site at CERN H2 beamline. The drawing 
shows various calibration sub-systems. 
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Fig. 7. H2(1995) data: ADC spectrum of 225 GeV/c muons in a single sampling 
scintillator layer. Based on the observed inefficiency of the layer of approximately 
20%, the average number of photoelectrons per minimum ionizing particle per layer 
is 1.6 pe/mip. 
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Fig. 8. H2(1995) data: average energy response of the tile/fiber calorimeter to pions 
and electrons as a function of B field, relative to response for B = Tesla field. Also 
shown on the plot is the relative response of scintillator to radioactive gamma ray 
calibration source, as a function of B field. Here the B field lines are perpendicular 
to the scintillator plates (endcap configuration). 
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Fig. 9. H2(1995) data: comparison of 300 GeV/c pion shower profiles for B=0 and 
B=3 Tesla magnetic fields. The B field lines are perpendicular to the scintillator 
plates (endcap configuration). The pion shower profiles are divided by the average 
muon response for each layer, which corrects for the overall scintillator brightening 
effect. 
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Fig. 10. H2(1996) data: comparison of 300 GeV/c pion shower profiles for B=0 and 
B=3 Tesla magnetic field. Here, the B field lines are parallel to the scintillator plates 
(barrel configuration). The pion shower profiles are divided by the average muon 
response for each layer, which corrects for the overall scintillator brightening effect. 
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Fig. 11. H2(1996) data: effect of a B field on the average energy response of the 
tile/fiber calorimeter to pions and electrons divided by the average muon response. 
Response ratio is normalized to 1 for B= Tesla. Also shown are the results from a 
GEANT simulation. Here, the B field lines are parallel to the scintillator plates (bar- 
rel configuration). The position of the scintillator package relative to the incident 
beam direction is the following: 1 mm plastic + 4 mm scintillator + 2 mm plastic 
(configuration A). The overall scintillator brightening B field effect is removed since 
the response of electrons and pions are divided by the average muon response in 
each layer. 
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Fig. 12. Orientation of the scintillator package with respect to the beam direction 
inside the absorber gap (configuration A) : 1 mm plastic + 4 mm scintillator + 2 mm 
plastic. The beam is incident from the left. In configuration B, the orientation of 
the package is reversed: 2 mm plastic + 4 mm scintillator + 1 mm plastic. 
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Fig. 13. Results of a GEANT simulation of the HCAL calorimeter response to 10 
GeV electrons (in a 3 Tesla magnetic field parallel to the scintillator plates) as a 
function of the air gap (d) between the scintillator package and the most upstream 
absorber plate for configuration A. 
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Fig. 14. H2(1996) data: effect of a B field on the average energy response of the 
tile/fiber calorimeter to electrons and pions. Response ratio is normalized to 1 for 
B= Tesla. Also shown are results from GEANT simulation. Here, the B field 
lines are parallel to the scintillator plates (barrel configuration). The position of 
the scintillator package relative to the incident beam direction is as follows: 2 mm 
plastic + 4 mm scintillator + 1 mm plastic (configuration B). The overall scintillator 
brightening B field effect is removed since the response for electrons and pions is 
divided by the average muon response. 
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Fig. 15. H2(1996) data: the ratio of the pion response of HCAL at 3 Tesla and 1.5 
Tesla to the pion response at Tesla versus beam energy. Here, the B field lines 
are parallel to the scintillator plates (barrel configuration). The data is shown for 
two configurations of scintillator package inside absorber gap: 1 mm plastic + 4 mm 
scintillator + 2 mm plastic (configuration A) and 2 mm plastic + 4 mm scintillator 
+ 1 mm plastic (configuration B). The lines correspond to the predictions using 
Groom (solid lines) and Wigmans (dashed lines) parametrizations of the fraction of 
the electromagnetic component in hadronic showers. 
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Fig. 16. H2(1995) data: energy deposited by 225 GeV/c muons in HCAL. The dashed 
line shows the energy reconstructed in HCAL for random triggers (pedestal events). 
The pedestal peak has a rms width of 80 MeV of equivalent hadron energy. 
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Fig. 17. H2(1995) data: scatter plots of Ehcal, the energy deposited in the 
hadron compartment, (vertical scale) versus Eecal, the energy deposited in the 
Pb-scintillator electromagnetic compartment, (horizontal scale). In addition to pi- 
ons, there is a contamination of electrons and muons in the beam. 
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Fig. 18. H2(1995) data: reconstructed energy distributions for " mip-in-EC AL" pi- 
ons. The minimum ionizing energy deposition in ECAL has been added to the energy 
sum. 
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Fig. 19. H2(1995) data: distributions of Etot, the sum of energies deposited in the 
hadronic and electromagnetic compartment for "full pion sample". 
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Fig. 20. H2(1995) data: linearity of the energy response for "mip-in-ECAL" pions 
and "full pion sample" . The statistical errors are smaller than the size of the symbols. 
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Fig. 21. H2(1995) data: fractional energy resolutions for "mip-in-ECAL" pions and 
"full pion sample". The statistical errors are smaller than the size of the symbols. 
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Fig. 22. H4(1995) data: response of the PDWO4 crystal ECAL calorimeter to 25, 
50, 100 and 150 GeV/c electrons. 
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Fig. 23. H4(1995) data: scatter plots of the energy in HCAL versus the energy in 
the PDWO4 crystal ECAL calorimeter for pions. 
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Fig. 24. H4(1995) data: the energy response of HCAL for "mip-in-ECAL" pions. 
The minimum ionizing energy deposition in ECAL has been added to the energy 
sum. 
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Fig. 25. H4(1995) data: the energy response of the combined crystal ECAL+HCAL 
calorimeter to "full pion sample". 
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Fig. 26. H4(1995) data: linearity of the energy response for "mip-in-ECAL" and 
"full pion sample". The statistical errors for some of the data points are smaller 
than the symbol size. 
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Fig. 27. H4(1995) data: fractional energy resolutions for "mip-in-ECAL" and "full 
pion sample" . The statistical errors for some of the data points are smaller than the 
symbol size. 
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Fig. 28. H4(1995) data: the response of HCAL to electrons. The ECAL crystals are 
moved out of the beamline for these data. The figure also shows the response of 
HCAL to pions before and after a correction for downstream longitudinal energy 
leakage is applied. A 30% systematic uncertainty is assumed on the longitudinal 
leakage correction. The dotted line represents the ratio of the average electron re- 
sponse to the pion response at 375 GeV/c. 
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Fig. 29. H4(1995) data: the pion/electron response ratio for HCAL (only) as a 
function of beam momentum. The pion response of HCAL has been corrected for 
longitudinal leakage. We assume a linear electron response for HCAL (only), and 
use E(ele)/E(7r)=1.15±0.01 at 375 GeV/c. 
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Fig. 30. H2(1996) data, B = Tesla data: the dependence of the E/p and the 
fractional rms energy resolution for 300 GeV/c pions as a function of parameter a, 
the weight assigned to HI (the first readout compartment of HCAL). 
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Fig. 31. H2(1996), B = Tesla data: the linearity of the response of HCAL to 
"mip-in-ECAL" pions and of the combined PbW04 crystal ECAL + HCAL system 
to "full pion sample" with a=1.4 and without the weighting of HI. The statistical 
errors are smaller than the size of the symbols. 
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Fig. 32. H2(1996), B = Tesla data: the fractional pion energy resolution of HCAL 
to "mip-in-ECAL" pions and of the combined response of the PDWO4 crystal 
ECAL+HCAL system to "full pion sample" with a=1.4 and without the weighting 
of HI. The statistical errors are smaller than the size of the symbols. 
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Fig. 33. H2(1996), B = Tesla data: average hadron shower profiles for 50, 100, 150 
and 300 GeV/c pions as a function of calorimeter absorber depth. The total Xint 
includes the contribution of ECAL. 
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Fig. 34. H2(1996) data, B = 3 Tesla: the fraction of 300 GeV/c pions with re- 
constructed energy less than 200 GeV (approximately 3 a below the mean, or 100 
GeV of missing energy) versus total absorber depth. The total Xint includes the 
contribution of ECAL. 
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Fig. 35. H2(1996), B = 3 Tesla data: comparison of the fractional pion energy 
resolution of the combined crystal ECAL+HCAL system, for various choices of 
absorber samplings in the inner HCAL:(a) 3 cm Cu sampling for first 8 layers, 
followed by 6 cm sampling, (b) 6 cm uniform sampling, and (c) 12 cm uniform 
sampling. 
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Fig. 36. Comparison of H2(1996), B = Tesla data with GEANT-GCALOR sim- 
ulations for the average longitudinal shower profile of 50 GeV/c pions. The electro- 
magnetic calorimeter has been taken out of the beamline. 
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Fig. 37. H2(1996), B = Tesla data: linearity of the energy response of the hadron 
calorimeter (HCAL only) for pions and comparison with the GEANT-GCALOR MC 
simulation. The electromagnetic calorimeter has been taken out of the beamline. 
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Fig. 38. H2(1996), B = Tesla data: the fractional energy resolution of the hadron 
calorimeter (HCAL only) for pions and comparison with the GEANT-GCALOR MC 
simulation. The electromagnetic calorimeter has been taken out of the beamline. 
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Fig. 39. H2(1996), B = Tesla data: the linearity of the energy response of the 
PbWC>4 crystal ECAL+HCAL combined system for the "full pion sample" and 
comparison with the GEANT-GCALOR MC simulation. 
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Fig. 40. H2(1996), B = Tesla data: the energy resolution of pions of the PbW0 4 
crystal ECAL+HCAL combined system for the "full pion sample" and comparison 
with the GEANT-GCALOR MC simulation. The MC simulation includes test beam 
detector effects (such as transverse leakage from prototype ECAL and ECAL elec- 
tronics noise). Also shown (triangle symbols) is the MC simulation of the energy 
resolution of pions in crystal HCAL+ECAL combined system without these test 
beam detector effects (which are not expected to be present in the CMS experi- 
ment). 
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Fig. 41. Comparison of GEANT simulation for the pion energy resolution of 
the crystal ECAL+HCAL combined system using various generators (GHEISHA, 
GCALOR, and GFLUKA-MICAP) of hadron showers. The test beam data is repre- 
sented best by the GCALOR simulation, which include all test beam detector effects 
and is labeled as "smeared". The GCALOR simulation labeled as "unsmeared" is 
expected to represent the experimental situation in the CMS experiment. 
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